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Abstract 

We give the lens equation for light deflections caused by point mass condensations in an otherwise 
spatially homogeneous and flat universe. We assume the signal from a distant source is deflected 
by a single condensation before it reaches the observer. We call this deflector an embedded lens 
because the deflecting mass is part of the mean density. The embedded lens equation differs from 
the conventional lens equation because the deflector mass is not simply an addition to the cosmic 
mean. We prescribe an iteration scheme to solve this new lens equation and use it to compare our 
results with standard linear lensing theory. We also compute analytic expressions for the lowest 
order corrections to image amplifications and distortions caused by incorporating the lensing mass 
into the mean. 

PACS numbers: 98.62.Sb 
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I. INTRODUCTION 

Conventional extragalactic gravitational lensing assumes that the Universe is homoge- 
neous and isotropic on scales significantly smaller than observer /source/deflector distances, 
i.e., that the cosmological principal applies at these distances. It also assumes that a lens- 
ing inhomogeneity such as a galaxy or cluster of galaxies is an addition to the homogeneous 
mean. What we investigate here is the extent to which errors are made because of this latter 
assumption. To assume a single galaxy is an addition to the mean might not seem irrational 
but to assume giant super clusters are is more suspect. In fact they are both contributing 
to the mean and hence do not act as infinite range deflectors. To understand why, one only 
has to surround a typical deflector by an imaginary sphere of radius r and note that the 
average mass density inside the sphere decreases as r increases until the density reaches the 
cosmological mean at some r—rj,- If this were not correct the cosmological principle would 
be in error. Beyond the gravitational boundary r^, the gravitational field has returned to 
the homogeneous mean and the lens ceases to produce any additional deflection of a passing 
light ray. In this paper we compute modifications to the lens equation caused by this finite 
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range. To make sure we properly account for the lensing gravity we use an exact solution to 
Einstein's equations. We assume the deflector is a simple point mass lens embedded in a flat 
Friedman-Lemaitre- Robertson- Walker (FLRW) universe, see Eq. (JTJ) , whose energy content 
includes pressureless dust (cold dark matter) and a cosmological constant A (fl TO + Q\ — 1). 
The mathematics of the embedding process is the same as embedding in the Swiss cheese 
cosmological models jl-4]. These models are the only known exact general relativistic (GR) 
solutions which embed spherical inhomogeneities into homogeneous background universes. 
The range r& above is given by the comoving radial boundary of the homogeneous sphere 
that has been replaced by the condensation. Beyond that boundary the gravity caused by 
a condensation and a homogeneous sphere are exactly the same. Schiicker j|| refers to this 
radius as the Schiicking radius. For a point mass lens the removed dust sphere of comoving 
radius Xb is replace by a Kottler condensation jg], i.e., Schwarzschild with a cosmological 
constant, see Eq. (|2J). In {7, 8] we derived analytical expressions for the bending angle a and 
the time delay AT of a photon that encounters such a condensation. Related work appeared 
in 9l-ll4j]. In this paper we derive the embedded lens equation and prescribe a scheme to 
iteratively solve it. 

The flat FLRW metric for the background cosmology can be written as 

ds 2 = -c 2 dT 2 + R(T) 2 [d X 2 + X 2 (d6 2 + sin 2 9d<p 2 )] , (1) 

and the embedded condensation is described by the Kottler or Schwarzschild-de Sitter metric 
[f| which can be written as 

ds 2 = - 1 {r)~ 2 c 2 dt 2 + 1 (r) 2 dr 2 + r 2 (d9 2 + sin 2 9 d(j) 2 ), (2) 



where 7 -1 (r) = Jl — f3 2 {r) and /3 2 (r) = r s /r + Ar 2 /3. The constants r s and A are the 
Schwarzschild radius (2Gm/c 2 ) of the condensed mass and the cosmological constant re- 
spectively. By matching the first fundamental forms at the Kottler-FLRW boundary, angles 
(9, 0) of Eqs. (T2J) and (DQ) are identified and the expanding Kottler radius r b of the void is 
related to the comoving FLRW radius Xb by 

r b = R(T) X b- (3) 

By matching the second fundamental forms the Schwarzschild radius r s of the Kottler con- 
densation is related to FLRW by 



n m %R oXb ) 3 , (4) 



where Hq is the Hubble constant and the cosmological constant A is constrained to be the 
same inside and outside of the Kottler hole. 

In Section II we give the lens equation valid for deflections caused by Kottler condensa- 
tions in the flat FLRW universe and numerically compare its predictions with conventional 
lensing theory for a source at redshift one and a deflector at redshift one half. In Section 
III we give analytic expressions for image magnifications and distortions for the embedded 
point mass lens (to lowest order only) and compare them with conventional lensing results. 

II. THE LENS EQUATION 

The Swiss cheese lensing geometry is shown in Fig. [TJ The deflected photon leaves a 
source S, enters a Kottler hole at point 1, exits at point 2 with a deflection angle a < 0, and 
then proceeds to the observer at 0. Point B is the intersection of the forward and backward 
extensions of respective FLRW rays SI and 20 drawn as if the Kottler hole were absent and 
the original ray was simply reflected at point B. Angles 9j and 9s are respectively the image 
and source positions relative to the observer-deflector optical axis OD. The rotation angle 
p measures the difference between the horizontal axis [with respect to which we measure the 
spherical polar angle 0, see Eqs. §2§ and ( ITS]) and Fig. 1] and the optical axis. A negative 
p is a clockwise rotation of the observer. The lens equation for a given deflector mass and 
background cosmology is simply the equation that gives 9j as a function of 9s for fixed 
comoving source-observer distance Xs and deflector-observer distance Xd, and fixed photon 
arrival time T . For non-embedded lenses, i.e., for conventional linear lensing theory, this 
relation is straightforward to obtain even for complicated lensing mass profiles, because the 
deflector is completely unrelated to the cosmology. For an embedded lens this is no longer 
the case. However, because of the azimuthal symmetry of the lensing geometry all photon 
orbit variables can be thought of as depending on a single independent variable. Choosing 
9s or the photon's minimum Kottler coordinate tq would be logical but not convenient. In 
what follows we have chosen to give all quantities as functions of where 7r — (f>i is the 
azimuthal angle of the photon at entry into the Kottler void (see Fig. 1, or Fig. 1 of {?]]). 
Because r ((fii) is a complicated function, r is retained in all expressions and only evaluated 
when needed. 

The embedded lens equation can be obtained by applying the law of sines to the triangle 
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FIG. 1. The comoving embedded lensing geometry. Points S, D and O represent respectively the 
source, deflector, and observer positions. The point B is a fictitious reflection point. Points 1 and 2 
denote the photon's entrance and exit from the Kottler void. The bending angle is a, 9j and 8s are 
respectively the image and source position angles at the observer measured relative to the optical 
axis OD, and A9 = 9j — 9s- A similar geometry appears in Fig. 1 of [Yj. The figure represents the 
9 = ir/2 plane (the plane containing the photon's orbit) of the spherical polar coordinates used in 
Eqs. (pQ) and (|2j). The <f> orientation is fixed by requiring the photon's point of closest approach to 
the Schwarzschild mass, vq, occur at <p = ir/2. 



SBO of Fig. 1 



sin (0s — 9 1 — a) 



Xbo 

Xs 



sin —a 



(5) 



The comoving distances \s and \d are often replaced by angular diameter distances D s and 
Da which are respectively functions of redshifts z s and za. The embedded lens equation can 
be compared to the standard linear lensing equation [15|, Q for flat FLRW 



9 s — 9 1 



■(-«) 



Xs Xd . 



(6) 



D s Xs 

where small angle approximations are made and the differences between distances from the 
observer to the deflector and to the reflection point B (xd and Xbo) are neglected. Since we 



are now computing the linear and non-linear corrections 



cannot make such simplifications as is done in 



;o the standard lensing theory, we 



17| and 



181 ] . To find the relation between 
these two distances we use the comoving triangle D20 and obtain 



X20 



cos 9i — cos(0i — £i + A0 — a 



Xd 



Xd, 



(7) 



where Xb, the Kottler void radius [see Eq. (Tj|], is assumed known. The angles £1, £2 = £1 +oc, 
0i, and A0 are exhibited in Fig. 1 and are the same as those used in [3, |8j where analytic 
expansions for them as explicit functions of r and 0i can be found. The angles £1 and £2 
are negative and give the respective slopes of the photon as it enters the Kottler hole at 
azimuthal angle ir — 0i and exits at angle 02 — A0 + 0i (see Fig. 1 of {?]]). The comoving 
distance XB2 can be obtained from trig identities applied to triangles 1B2 and 1D2 of Fig. 1 



XB2 



o sin(-A0/2 + £i) 
-2 cos 



A0\ 
1 + — I Xb- 



sin a \ 2 j 

Combining this with Eq. ((7j) we obtain the relation of xbo to Xd, 



(8) 



Xbo = X20 + XB2 
= I cos Qi — 



cos(0i — £1 + A0 — a) 
A0 S 



2sin(-A0/2 + £1) 

H cos 

sin a 

g{9 I ,^ 1 ,A<f>,a)xd- 



+ 



Xb_ 

Xd 



'Xd 



(9) 



The new lens equation ([5]) becomes 

s = e I + a + sin" 1 



Xd 

Xs 



g(9i,£i, A0,a) sin(-a^ 



(10) 



The task at hand is to evaluate all variables on the right hand side of Eq. ( flOl) as functions 
of a common variable e.g., 0i. Once accomplished, 6s{<pi) and #r(0i) can be tabulated to 
give the desired image position as a function of the source position, 9i(9s). The image angle 
9 1 can be determined from knowledge of £1, A0, and a by applying the law of sines to the 
triangle D20 



sin 6 1 / = sin(0i — £1 + A0 

a, 



a) 



Xb 



(11) 

is given by Eq. (13) of jg], 



Xd 

The bending angle a is given by Eq. (32) of (7J, A0 = 2 — <p\ is given Dy 
and the photon's slope angle £1 results from evaluating Eqs. (16)-(19) of [7| at the photon's 
entry point into the Kottler void (to fourth order) 

~ 7TI ~ ~ \ Tfl 

ii = —fi\ sin 0i H cos 0i (2 + sin 2 X ) - -fix — (6 

r 3 r 

~ 1 ~ 1 "IT? 

-3 sin 2 0i-2 sin 0i) - — ^Ar^ sin X - -— 

18 4 Tq 

x 15(01 — — ) + cos 0i (8 — 15 sin 0i + 4 sin 2 0i 



+14 sin 3 0i + 4 sin 5 0i) 



+ 0(5). 



(12) 
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The rotation angle p can be computed from the photon's exiting slope £1 + a and the image 
position 9j using 

p = £1 + 01 + 0! <0. (13) 

The expansion speed (5 = v/cof the void boundary relative to stationary Kottler observers 
is defined in Eq. (2) and when evaluated at the photon's entry point is called fti (sec Fig. 1 
of ~\). Keeping terms to 4 th order is necessary in order to correct point mass time delays 
for embedding. 

In the expressions for £i, A0, and a, approximation orders have been counted as follows: 
j3i is 1 st order, r s /ro and Ktq are both 2 nd . All terms are made of sums and/or products of 
these. The expansion speed f3\ depends on 0i and r through its dependence on t\ (which 
is given by the symmetric null geodesies of the Kottler metric Eq. fl2])) 

r o L , r s / , . ~ 2 \ / r s \ - 



rx=^ T <1 + 7T- 1 + sin^--^ - x 
sin 0i { 2r \ sin 0i / V2r / 

IT 1 2 7 4 15 



- - sin 0i - , - + — ( 7T - 20i 1 cot 0i 

4 4 sin 0i o v 7 

+0(6)J. (14) 

The above expansion is valid only when sin0i ^> r s /r . All quantities on the RHS of the 
embedded lens equation ffTUj) can now be evaluated as functions of 0i and r . These two 
variables fix the photon's symmetric orbit (symmetric about = tv/2) while in the Kottler 
hole. They are independent unless the photon is additionally constrained by originating at 
a specific cosmic source or arriving at a specific observer. To eliminate one of these two 
variables an additional relation between them such as a cosmic timing constraint must be 
used. For the photon which started at a fixed Xs to have reached the observer at time T 
after entering the Kottler void at 0i and passing with minimum impact ro, it must have 
impacted the Kottler void at a specific time T\ or equivalently at a specific redshift z\ 
(1 + z\ = Rq/ R{T\)). Knowledge of z\ allows us to independently determine T\ from the 
embedding equations (EJ) and (jlj) i.e., by using 

ri= iT^(jiwJ ■ (15) 

Because Z\ is not assumed known we compute Z\ — z^, the difference in entry redshift and 
the (assumed known) deflector redshift, using techniques similar to those developed in [7|, |8|. 



The result up to fourth order is 



where 



and 



In the above 



z i -z 1 = (l + Zl ) A;? lst (zi,0i) + A^ nd (z 1? r 



> VI, 



+Az 3rd {z 1 ,r ,<f> 1 ) + Az 4th (*i,r o ,0i) 



4th, 



Az lst = COS 0i, 



Az 



2nd 



1 m 



Ar , 1 a Xb . 2 

— — + sm (f>i + - — sin pi 

3 2 Xd 2r 

—7 sin 2 0i 



(3 



Ar Xh rim 

— cos 0i + -ft — 

3 3 r 



Az 3rd = i/3iAr o 2 cos0i-^ 

3 

/ 7 r 72i <M ~ 

cos 0i 7 + 26 sin 0i + 12 log tan — j sin0i 



x 



7m Xb 7-37 

cos 0i sin 0i, 

2r Xd 



„4th 



A: ;ih = bM-a - 2sin0 2 ) + Ui^-U + 9sin0 2 
6 Xd 2 Xdr K 

1 / \ 3 



1 • / Xb \ . 4 , 

— ) sin 0i + - — 



3 m 
!r 



— 18 sin 01) sin 0i + -j3i 

8 \Xd, 

Ar 2 , esc 0i (61 + 24 sin 0i 

36 r v 



x I — I sin 5 0! - — — 
\XdJ 36 r 

+ 124 sin 2 - 227 sin 4 + 48 cos 0i log tan y 

+ —^ (36 - 18 sin 0i -431 sin 0? 
12 rg v 

+ 42 sin 0? - 188 sin 4 + 595 sin 0? 

or 01 \ 

— 240 cos 0i sin 0i log tan —J . 



Xb = 



1 /r. c 2 \ V3 1 
1 + z d \Q m H% J D d 

larger 



than 1 st order 
;ing fc 



(16) 



(17) 



(18) 



(19) 



(20) 



four variables Z\, r , ri, 
and ri as 



(21) 
and 



Xd L^r z d \ 

is taken as an additional small parameter no larger man 1 

Equations (TTBl) . ( TT4"1) and (fl5l) are three equations relating 
0i. They can be solve iteratively (four iterations) giving z±, tq, and r\ as functions of 0i. 
For an example, to obtain z\ correct to the first order in we use Eqs. ffTB"]) and ([T 



Z 1 = Z d —(l + Z d )Az S 0i) = Z d + (1 + Z d )/3(z d ) COS( 



'i- 



(22) 
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FIG. 2. The embedded point mass lens versus the Schwarzschild lens. The deflector /source redshifts 
are respectively = 0.5, z s = 1.0; the cosmological parameters are fL m = 0.3, JIa = 0.7, and 
Ho = 70kms~ 1 Mpc^ 1 ; and the deflector mass is m = 10 15 M@, The abscissa 9s is the source 
angle measured in units of the classical Einstein ring angle 9e and the dashed/dotted lines are for 
primary /secondary images. The bifurcating blue curves are above the corresponding bifurcating 
red curves. The green bifurcating pair of curves in the right panel are between the upper blue pair 
and lower red pair. The solid curve in the left panel measures the relative correction of the angle 
between the primary and secondary images. 

this can be inserted into Eq. (1151) to obtain r\ correct to first order in This r\ is then 
inserted into Eq. ( jb£l) (only the lowest order is needed here, i.e.,r = risin0i) to obtain 
r correct to first order. For the next iteration, we include Eq. (|T8|) and the r s /r term in 
Eq. (fT^j) . and so on. With Zi(za, <f>i), r (zd, and r 1 (2 d ,0 1 ) in hand, we can compute 
£i, A(j), and a in terms of <pi and finally solve the embedded gravitational lensing equation 
( fTOl) for 6s(4>\) which can be tabulated to give 6s(0i) for a given image. 

In Figs.[2]and[3]we have solved the embedded point mass Swiss cheese lens equation (JTUJ) 
and compared the results with those of the conventional Schwarzschild point mass lensing 
theory. We chose deflector /source redshift respectively Z& = 0.5, z s = 1.0, cosmological 
parameters D, m = 0.3, = 0.7, and H = 70kms~ 1 Mpc -1 . In Fig.El we chose a deflector 
mass m = 10 15 M (a rich cluster). For each source angle 9s, we solved Eq. ( jTOj) using the 
iteration scheme described above obtaining <pi, Z\, tq, n, 9i, etc., for both the primary and 
secondary images. The conventional Schwarzschild results are given by Eq. ([6]). The impact 
parameter in conventional lensing is simply taken as r ( Sch ) = 9i(s c h)D d . The dashed/dotted 
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FIG. 3. The embedded point mass lens versus the Schwarzschild lens. Same as Fig.[2l except that 
the deflector mass m = 10 12 Mq. 

curves are for primary /secondary images, and the solid curve is the correction to the angle 
between image pairs, i.e.,#n — 9n- In the left panel, we compute the relative correction in 
the image position, i.e.,S9j/9i(Sck) (blue-upper bifurcating pair of curves), and the relative 
correction of the impact parameter r , i.e., 5r /r (s c h), where 5r Q = r — r (s c h) (red-lower 
bifurcating pair of curves). 

In the right panel, we compute the net correction of the bending angle a (central pair 
of green curves), the effect of the linear correction alone, i.e., cos 3 0i — 1 (lower pair of red 
curves), and the contribution of the cosmological constant A (upper pair of blue curves). 
Figure |3] is the same as Fig.[2] except that it is for m = 1O 12 M (a typical large galaxy). For 
m = 10 15 M , corrections in the image angle 9j can be as large as 0.2%, and corrections in 
the bending angle a can be as large as —0.8%. For m = 10 12 M , corrections in the image 
angle 9j can be as large as 0.01%, and corrections in the bending angle a can be as large as 
-0.18%. 



III. IMAGE MAGNIFICATION AND ELLIPTICITY 

In this section we include only the lowest order correction to the standard lensing equa- 
tion caused by the finite range of the embedded point mass Swiss cheese lens. Sereno [19 1 
computes alterations in the magnification but only within the Kottler void. We assume 
sin#/ <C 1, sin 6*5 <C 1, and that the Kottler hole is much smaller than the observer-deflector 
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distance, i.e., g(<pi) -1 < 1, see Eq. flO]). From Eq. (jSJ) we obtain 

0s-0i = -^{-a), (23) 

which is the same as the standard lens equation (jHJ) except that the bending angle to the 
lowest order now contains a cos 3 0i factor caused by the finite range of the deflector 

a = -2— cos 3 0i, (24) 
r 

see Eq. (32) of p|. Equation (jSJ) is the form assumed correct by [3] but with a different 
expression for the deflection angle a. 

To lowest order the minimum Kottler impact is 

r = D d e I + O((3 1 ), (25) 

[see Eqs. (j3J), (fill and (fl4j) ] and the embedded lens equation to lowest order becomes 

o s -ei = - d -jr cos 3 0i. (26) 

Pi 

The angle #e is the familiar Einstein ring radius 



(27) 



and from Eq. ( ffTT) 6 1 / is related to 0i by 



This gives us a modified Einstein ring radius (to lowest order) 



sin 0i = + 0{p x ). (28) 



^ = V 2 ^S (co ^ l)3/2 ' (29) 

(see 20] for modifications in the Einstein ring within the Kottler void). The two images for 
the standard point mass lens are easily found at 

3± 1 



0f=2(0.±V*2+40!j, (30) 

however, to find the corresponding image positions for the embedded lens you must solve 
Eq. flU with 6 E replaced by 6' E . 
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The amplification and shear for the embedded lens can be found by a familiar 15] rescaling 
9s — > 0s/0e = y, di — ► 6i/0e = Equation f l26|) simplifies to 



where 



y = x 



sin (pi 



COS 3 0i 
X 2 



The 2-d Jacobian A = <9y/<9x is found to be 16j 



A 



COS </>i 



1 
1 



cos 0i (2 + sin^ 



x 1 



2 

X\X<2 



X±X 2 X2 



which has two eigenvalues 



(31) 



(32) 



(33) 



a± = 1 + cos 0i (1 + 2 sin 2 0i)— ~ 



x 



2 ■ 



, ~ 1 

a 2 = 1 - cos 0i^t. 

or 



Writing 



A 



(34) 



(35) 



' 1 - K - 71 -72 \ 

V -72 1 - « + 7i / 
as is commonly done in standard gravitational lensing theory, we immediately obtain a 

negative surface mass density 



and two shear components 



3 • 21 I 1 
k = --sin 0i cos 0i — , 

2 a;" 1 



iy 2 ™ 2 

Tl = -cos 0i (2 + sin 0i ) 2 ^ 4 ■ 
7 2 = -cos 0i (2 + sin 0i) — — , 

37 



with total shear 



7 = V7i +72 = cos0i(2 + sin 2 0i)- o ,. 



1 

2x~ 2 



The amplification /i for an image is given by 

/i _1 (x) = det A — (1 — k) 2 — 7 2 = ctici2 



T o T 1 /IT/ n "* . 1 

l + 3cos0isin 0i— — cos 0i(l + 2sin 0i; 



(36) 



(37) 



(38) 
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FIG. 4. Linear corrections to Schwarzschild lensing caused by the finite range of embedding — the 
magnification ratio Hi/ ^2 and the ellipticity e are plotted as a function of source position. The 
cosmological parameters and redshifts are same as in Figs. 2 and 3. 

The image of a circular source (eccentricity e = 0) will be an ellipse of eccentricity 



a-2 y(2x 2 + 3sin <j>\ cos 0i)(2 + sin 0i)cos0i 
^ a\~ x 2 + cos0i(l + 2sin 2 1 ) 

The standard lensing results are obtained by putting cos0i = 1 and sin^i = in the 
above. Deviations from standard image amplification \i and the image ellipticity e caused by 
embedding are shown in Fig.Hl The left panel is for a deflector mass m = 10 15 M and the 
right is for m = 10 12 M & . In each plot, the red solid and the (identical to accuracy shown) 
black dotted (upper) curves show the corrections in ellipticity, i.e., 5e/e for the primary and 
secondary images. The solid blue (lower) curve is the relative correction in the magnification 
ratio, i.e., 5(fii/ ^2) / (Hi/ ^2)- F° r the m = 1O 15 M case, the correction in ellipticity can be as 
large as 0.03%, and the correction in magnification ratio can be as large as —0.17%. For the 
m = 1O 12 M case, the correction in ellipticity can be as large as 0.004%, and the correction 
in magnification ratio can be as large as —0.019%. 



IV. CONCLUSIONS 



We have given a lens equation (jHJ) valid for use on highly concentrated lenses (point 
masses) which are embedded into the otherwise spatially homogeneous and flat background 
FLRW cosmology. We have also given the additional equations necessary to iteratively solve 
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this embedded lens equation and have outlined a procedure for doing so. As an example we 
have looked at differences in strong lensing predictions made by this new theory as compared 
to the conventional theory. We used a large galaxy size lens (m = 1O 12 M ) and a rich cluster 
size lens (m = 10 15 M & ) and found, as was suggested before in 7, Q], that predictions for 
strong lensing effects made by embedded lens theory differs by less than 1% from predictions 
made by the conventional theory. In Section II we looked at image angle differences and in 
Section III we looked at lowest order analytic expressions for image differences. We expect 
more significant effects to occur for weak lensing applications where impact distances are 
much larger and where shielding effects (cos 3 <j>i) are more significant. 

Work on this project was partially supported by NSF grant AST-0707704 and US DOE 
Grant DE-FG02- 07ER41517 and B. Chen wishes to thank the University of Oklahoma 
Foundation. 
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